Simple model for the spherically- and system-averaged pair density: 

Results for two-electron atoms 
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As shown by Overhauser and others, accurate pair densities for the uniform electron gas may 
be found by solving a two-electron scattering problem with an effective screened electron-electron 
repulsion. In this work we explore the extension of this approach to nonuniform systems, and we 
discuss its potential for density functional theory. For the spherically- and system-averaged pair 
density of two-electron atoms we obtain very accurate short-range properties, including, for nuclear 
charge Z > 2, "on-top" values (zero electron-electron distance) essentially indistinguishable from 
those coming from precise variational wavefunctions. By means of a nonlinear adiabatic connection 
that separates long- and short-range effects, we also obtain Kohn-Sham correlation energies whose 
error is less than 4 mHartree, again for Z > 2, and short-range-only correlation energies whose 
accuracy is one order of magnitude better. 



I. INTRODUCTION AND SUMMARY OF 
RESULTS 

Density Functional Theory (DFT) 0, S is nowa- 
days the most widely used method for electronic struc- 
ture calculations, in both condensed matter physics and 
quantum chemistry, thanks to a combination of low com- 
putational cost and reasonable accuracy. 

In the application of this theory within the Kohn-Sham 
(KS) formalism one deals with a model system (the 
KS system) of N noninteracting electrons in a local po- 
tential vks ( r ) that forces them to yield the same density 
n(r) of the physical system. The energy of the physi- 
cal system is then obtained from that of the KS system 
via a functional of the density, whose only term not ex- 
plicitly known is the exchange-correlation energy E xc [n]. 
Correspondingly, in the local potential vks (r) there is an 
unknown term, v xc (r) = SE xc [n]/5n(r). 

The success of KS DFT is mostly due to the fact that 
even simple physical approximations of E xc [n], like the 
local density approximation (LDA) £|, already give ac- 
ceptable results for many purposes. This spurred fun- 
damental research in the field, and led to a wealth of 
more and more sophisticated exchange-correlation func- 
tional 0, S 0] i an d to the development of different ap- 
proaches to DFT 0, ■ 

Recently, in the search for accurate E xc [n], the fo- 
cus of a large part of the DFT community has shifted 
from seeking explicit functionals of the density like the 
generalized gradient approximations (GGA) [8j, to im- 
plicit functionals, tipically using the Kohn-Sham orbital 
kinetic energy density (3 or the Kohn-Sham orbitals 
(see, e.g., UTEJEI)- The so-called "third generation" 
of exchange-correlation functionals is based on the ex- 
act exchange of the noninteracting (KS) system, sim- 
ply obtained by putting in the formal expression for the 
Hartee-Fock exchange the Kohn-Sham orbitals <pi a {r). 
Such expression corresponds to an implicit functional 
of the density, E x [n] = E x [{(pi a [n]}]. The local poten- 
tial v x (r) — SE x [n]/Sn(r) that generates the orbitals 
(fia- (r) [n] can be obtained via the optimized effective po- 



tential method (OEP) Qj. 

In this broad context, sketchily summarized here, we 
propose a simplified method to build the "bridge" be- 
tween the physical and the KS system, or, more gener- 
ally, with a reference model system of partially interact- 
ing electrons. We focus on a quantity which is known to 
play a crucial role in DFT and has an intuitive physical 
meaning, the spherically and system-averaged electronic 
pair density f(r 12) (also known in chemistry as spherical 
average of the intracule density, see e.g. |13L ll4L I15L Il6fl . 
and especially PJlllS)- Given the spin-resolved diagonal 
of the two-body reduced density matrix, 



7^l 2 (ri,r 2 ) = J \^{ri^i,...,r N a N )\ 2 dr 3 ...dr 



N, 



(1) 



we define the spin-summed pair density ri2(ri,r2), 

n 2 (r 1) r 2 ) = ^^V 7 ( 2 ) (ri ,r 2 ), (2) 



and we integrate it over all variables but ri 2 = |r 2 — 1*1 1 
by switching, e.g., to ccnter-of-mass coordinates, R = 
|(ri + r 2 ), r i2 = r 2 - ri, 
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The function /(ri 2 ) times the volume element 47rr^ 2 dr 12 
is proportional to the probability density for the particle- 
particle distance in a system of N electrons in the state 
ty, and is normalized to the number of electron pairs, 
N(N — l)/2 . This quantity fully determines the ex- 
pectation value of the electronic Coulomb repulsion (in 
Hartree atomic units used throughout), 
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(4) 



and is a measurable quantity, being essentially the 
Fourier transform of the electronic static structure fac- 
tor ^^|. By construction, the one-electron density n(r) is 
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FIG. 1: Spherically and system-averaged pair densities for two-electron atoms: 'exact' results | 4C| are compared with the values 
obtained for the Kohn-Sham system and with the present approach, which is designed to get realistic /(V12) starting from the 
Kohn-Sham ones. 



the same in the KS and in the physical system, whereas 
/(712) will be different in the two cases, as shown, e.g., in 
Fig. n for some two-electron atoms. In the physical sys- 
tem /(712) has a much lower "on-top" value /(ri2 = 0) 
than in the KS system, and it has a cusp |2(j , as expected 
from the fact that the electrons repel each other via the 
Coulomb interaction. Roughly speaking, in the classic 
DFT approach to correlation, the difference in energy 
arising when we evaluate the r.h.s. of Eq. J3J with the 
two f{r 12), the physical and the KS, is what one tries to 
describe with a universal functional of the density |2l| . 
Here we follow a different approach: we try to build re- 
alistic f{r\2) from a set of simple radial equations, to be 
solved for each system, and eventually coupled to a DFT 
calculation. 



formalism (Sec. [nj, we give a physically- motivated pre- 
scription for the effective e-e potential fSec. IHIfl . and we 
test it on the simple but not trivial case of two-electron 
atoms (Sec. llVfl . The prescription for the effective e-e 
potential used here is not very sophisticated. Improve- 
ments along the lines of what has been done for the uni- 
form electron gas [U 0,113 will be the subject of future 
work. Yet, even at this simple first stage of the theory 
we already obtain rather accurate results, especially for 
the short-range part of f[rv2) (see Fig. ^ an d Table QJ. 
In Sec. we show that with the present approach we can 
also recover the difference in kinetic energy between the 
physical and the KS system. Finally, Sec. is devoted 
to conclusions, perspectives and open questions. 



Our approach is inspired by the seminal work of Over- 
hauser [23] and its subsequent extension j2^| , in which the 
function f{r 12) for the uniform electron gas is obtained 
from a set of geminals, solutions of a radial Schrodinger 
equation with an effective electron-electron (e-e) poten- 
tial. Simple approximations for such effective e-e po- 
tent ial g ive indeed accurate results at all relevant densi- 
ties [23, Hj, |2fj ■ Here we try to generalize this approach 
to systems of nonuniform density to get accurate f(r 12). 
The main goal of the present work is understanding 
whether the method is promising, and whether it is worth 
developing and refining it. To this purpose, we define the 



II. FORMALISM 

In addition to the work on the "Overhauser model" jU 
I23I I2H , the approach described here takes advantage of 
inspiring papers on the possibility of constructing a pair- 
density functional theory [26l I27I and a local- 
density-of-states functional theory |3C|. 

Our starting point is a constrained search over ^N(N— 
1) "effective" orthonormal geminals ipi(r 12) that mini- 
mize the electron-electron relative kinetic energy T±2 = 
— V^ 12 (the reduced mass for the relative motion is 1/2) 
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and yield the exact /, J2i \^i( r i2)\ 2 = f(r 12 



(5) 



thus leading to a set of radial equations formally similar 
to the KS ones, 

[-Vr 12 + v c s(r 12 )}ipi(r 12 ) = Ci ipi(r 12 ) (6) 
E l^(^i 2 )| 2 = /(ri 2 ). (7) 

These equations imply that an expansion in spherical 
harmonics of f{r\ 2 ) has been done, so that the kinetic 
energy operator also contains the usual £(£ + l)/r 2 term. 
To fully define these equations we need a rule for the oc- 
cupancy of the effective geminals. In analogy with what 
has been done for the uniform electron gas [2H, |24| , we 
can assign spin degeneracy 1 to even-angular-momentum 
states (singlet) and spin degeneracy 3 to odd-angular- 
momentum states (triplet), up to N(N — l)/2 occupied 
states. More generally, for open shell systems it could 
be better to develop the formalism for the spin-resolved 
quantities, starting from Eq. 0. This will be investi- 
gated in future work. 

The effective electron-electron potential v c g(ri 2 ) of 
Eq. is the Lagrange parameter for f(ri 2 ), and is a 
functional of / itself and of the electron-nucleus external 
potential V ne (or, equivalently, of the density n(r)). To 
see this, we can rewrite our Eqs. ©-0 m terms of a 
minimization of the total ene rgy in two steps, using the 
constrained search formalism |3lll32j| for the ground state 
energy E = min*(#|T + V ee + V ne \^), 



E 



mm 

/ 



V 2 12 |^)}.(8) 



mT + v ne m 



Defining the kinetic and external-potential functional as 

^ke[/; Vne] = 
wmMT + Vnel*) - min V(^l - V 2 12 |^), (9) 



we can rewrite 



E = min 

/ 

+F+ 



{ , n \ in , E^i - v 'i, i^) + / — dr i2 

:e[/;K £ ]}. (10) 



Searching this minimum by directly varying the ipi (with 
given, fixed, V^ e ) leads to Eqs. ©-0 with the identifi- 
cation 

, v 1 , SF KE [f;V ne ] 
VeS{ri 2 ) = — + — JJJ 7 — • (11 



Thus, in principle we could recover the whole ground- 
state energy via the (unknown) system-dependent func- 
tional Fke[/; V ne ]. In practice, it seems much more fea- 
sible to combine Eqs. ©-(0 with a DFT calculation, 
that yields the complementary information (the density, 
and thus (*|V^ e |*))- The steps of Eqs. can be 

repeated for arbitrary electron-electron interaction and 
external one-body potential. In particular, we can set 
V e e = ^Vee and V ne — V x , where V x is an external po- 
tential that keeps the density equal to the one of the 
physical system. One could thus obtain f x at each cou- 
pling strength A between and 1 from Eqs. @-0 with 
a suitable v x s . The correlation energy of KS theory is 
then simply given by 0, 0] 

EM-f^S^ J^-'^ - (12) 

JO J r 12 

Alternatively, this procedure (usually called adiabatic 
connection |34|) can be performed along a nonlinear path, 
e.g., by setting [|| |3^, |3(| |3l} v x e = erf(Ar)/r, where 
erf (a;) is the error function (see Sec.[VJ. Eventually, the 
two sets of equations, KS and (©-0) plus 112(1. could 
be solved together self-consistently. This last issue is dis- 
cussed in Sec. lVIl Notice that if we combine Eqs. ©, 
and l|12f> with a DFT calculation, we only need to approx- 
imate the potential v x s (ri 2 ) and not the whole functional 
i*KE since the remaining information is provided by DFT. 

It is also worth to stress at this point that there is no 
wavefunction behind our Eqs. 0-0: the effective gem- 
inals ipi are defined via Eq. 0, and by specifying their 
occupancy (e.g., triplet and singlet). A bosonic version 
of the theory, in which only one geminal (proportional 
to yj f(r 12)) is occupied can also be considered [2^, |38| . 
In this work we only focus on two-electron systems for 
which the two choices are equivalent. A careful com- 
parison of performances of the "fermion-like" and of the 
"boson-like" occupancy in the uniform electron gas is the 
subject of current investigations 39]. 

As for KS DFT, the formalism just described can be 
useful only if simple approximations for v e ff(ri 2 ) yield ac- 
curate results. This is what we start to check in the rest 
of this paper. First, we construct a physically-motivated 
Ueff for two-electron atoms for the fully-interacting sys - 
tem, and we compare our results with "exact" ones |40j . 
Then, we generalize our construction to build v c g along 
the adiabatic connection, and we calculate the KS corre- 
lation energy. 



III. EFFECTIVE ELECTRON-ELECTRON 
POTENTIAL: THE OVERHAUSER MODEL 

For the interacting electron gas of uniform density n, 
Overhauser (2^ proposed a simple and reasonable effec- 
tive potential f e ff (^12): he took the sphere of volume n~ l 
around a given electron as the boundary within which 
the other electrons are excluded, due to exchange and 
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correlation effects. In the standard uniform-electron-gas 
model, a rigid positively-charged background maintains 
the electrical neutrality. Thus the exclusion region (or 
"hole") around a given electron, modeled with a sphere 
of radius r s = (47m/3) -1 / 3 , uncovers the background of 
positive charge, leading to an effective screened Coulomb 
potential with screening length r s , 



^(r ia ;r.) 



1 



W\<r. 



|r-ri2| 



dr, 



(13) 



equal to 



i 

''12 



3 
2r. 







n 2 > r a . (14) 



Equations ©-(01 combined with the Overhauser effec- 
tive potential of Eq. 1|14|) gave extremely accurate results 
for the short-range part (ri2 < r s ) of the function f(rW) 
in the uniform electron gas at all relevant densities |23| . 
A more sophisticated effective potential, based on a self- 
consistent Hartree approximation, extended such accu- 
racy to the long-range part of /(r , 12 ) at metallic den- 
sities |24j. Other approximate u B ff (fi?) for the uniform 
electron gas have also been proposed [25j , and exact prop- 
erties have been derived |41j. 

To produce realistic f(r 12) for nonuniform systems 
from Eqs. ©-0, here we generalize the original idea 
of Overhauser [22, Hj| to two-electron atoms, and show 
that it gives rather accurate results, especially for the 
short-range part of f(ri2). We start from the effective po- 
tential Vtff (ri 2 ) that generates /Ks( r i2)) the spherically- 
and system averaged pair density of the Kohn-Sham sys- 
tem. In the special case of a spin-compensated two- 
electron system, the KS wavefunction is simply equal to 
\\Jn{Yi)^Jn{r2)- Because, at this first stage, we are in- 
terested in testing our method as a "bridge" between 
the KS and the real system, here we use the "exact" 
Kohn-Sham system. We thus take accurate one-electron 
densities |40j . and construct /Ks( r i2), 



/Ks(fl2) 



n (R- V -y) n (it- 



2 J 4tt 



(15) 

and the corresponding "exact" potential - (ri^), that 
can be calculated by inverting Eqs. ©-0, 



,(0) _ 




const. 



(16) 



For systems with more than two electrons, the poten- 
tial v !S could be calculated, e.g., with the methods of 
Refs. |42ll43| . In practice, it would be much more efficient 

„(°) 



to build approximations also for v c ^ (see Sec. IVI|I . Ex- 
amples of functions /ks f° r nuclear charges Z = 1, 2, 3, 4 
are given in Fig. ^ they have a maximum at r 12 = 0, 
as expected in a system of two non-interacting electrons 
with antiparallel spins in a confining one-body external 
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TABLE I: Our results for the function /(n2) for two-electron 
atoms (first line for each property) compared with the corre- 
sponding "exact" quantities |4fij |. In the first line of the table 
we report the average r s as defined by Eqs. 1171 and I19H . For 
the "on-top" value /(0) we also show the LDA result (with 
/(0) for the uniform electron gas from Ref. |23|). All values 
are in Hartree atomic units. 



potential. When the interaction is turned on, the aver- 
age distance between the two electrons increases, with 
the constraint that n(r) is kept fixed. We can thus imag- 
ine that, with respect to the Kohn-Sham system, in the 
physical system the Coulomb repulsion between the elec- 
trons creates, on average, a screening "hole" around the 
reference electron of volume (7T) -1 , where n is an average 
density (i.e., n(r) integrated over the wavefunction), 



■ft J drn{r) 2 . 



(17) 



An approximate w ff(' , i2) could thus be simply con- 
structed as 



v cS (r 12 ) » v^(r l2 )+v^(r u ;r s ) 
with an average f s in of Eq. (|14fl , 



/47T 

v 3 



-1/3 



(18) 



(19) 



The Overhauser-like potential V°g(ri2', r s ) is thus a cor- 
relation potential to be added to the one that generates 
/ks- It describes the correlation between pairs of elec- 
trons due to Coulomb interaction, and keeps the informa- 
tion on the one-electron density in an approximate way, 
via the average n of Eq. jlTjl. Of course, for more com- 
plicated systems we expect to need a more sophisticated 
construction for r% . 



IV. RESULTS 

We have inserted the potential of Eq. (|18|l into Eqs. ©- 
0, and solved them for several two-electron atoms. Our 
results are shown in Fig. ^ an d summarized in Table [I] 
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We see that the simple effective potential of Eq. ff 1 81) gives 
already reasonable results for Z = 1 and 2, and that the 
accuracy of the results increases with Z (as the system 
becomes less and less correlated). The "on-top" value 
/(0) is essentially exact for Z > 2, and is much better 
than the LDA estimate (normally regarded as accurate) 
for all Z . This feature is appealin g, si nce the on-top value 
plays an important role in DFT |44|, and accurate /(0) 
are not easy to obtain from ab initio methods (see, e.g., 
Ref. m| and references therein). The term l/r 12 in the 
effective potential ensures that the calculated f(r 12) sat- 
isfies the exact cusp condition /'(0) = /(0). Table|l]also 
shows that the position r 1 ^* and the height f(r™ 2 x ) of 
the maximum of / is very well predicted by the present 
approach. The presence of this maximum is essentially 
due to the combined effect of the Coulomb repulsion be- 
tween the electrons and the confining external potential. 

In Fig. |21 we consider He and Ne 8+ , and we compare 
the correlated part of our /, f c = f — /kSi with the 
"exact" result |4(| and with the corresponding quantity 
calculated within LDA, i.e., 

/c LDA (n 2 ) = \J n(r) 2 g c (r 12] n(r)) dr, (20) 

where g c is the pair-correlation function of the uni- 
form electron gas at full coupling strength, taken from 
Ref. 0| . (For an extended system of uniform density n, 
we have g c = 2f c /nN.) Figure shows the same quanti- 
ties multiplied by 47rri 2 , i.e. the integrand of Eq. Q for 
the correlation part of (V ee ): the area under each curve 
gives (V ee ) — (V ee ) /fs- In the last line of Tabled we report 
quantitative results for (V e e) — (V ee )KS- This quantity is 
less accurate than the short-range properties, but it is 
still encouraging. Moreover, it saturates for large Z as in 
the exact case. 



V. ADIABATIC CONNECTION AND 
CORRELATION ENERGY 

For the calculation of the energy of the physical sys- 
tem, in addition to V c [n] — (V ee ) — (V ee )KS, one needs to 
know the kinetic- energy difference, T c [n] = (T) — (T)ks, 
that can be obtained via the adiabatic connection for- 
malism 0, l33l l34j . By varying a parameter A, the in- 
teraction v x e (ri 2 ) between the electrons is switched on 
continuously from zero to l/ri 2 , while the density is 
kept fixed by an external one-body potential V x . If 
v x =° — and v x = a — l/ri 2 , the KS correlation energy 
E c [n] = T c [n] + V c [n] is given by 0,111 

E c [n] = £ dX J™ dr 12 4tt r\ 2 / c A (r 12 ) dv *j r ^ , (21) 

where / C A = / A - / KS - 

Usually, the adiabatic connection is performed along a 
linear "path" 0,H3, by setting v x e = A/r 12 , which leads 
to Eq. I|12|) . If one is able to compute the exact f x , the 
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FIG. 2: The correlated part of the spherically- and system- 
averaged pair density, fc.(j\i) = f{r 12) — /ks(fi2)- Our results 
for He and Ne s+ are compared with the "exact" ones and with 
the LDA result (the hole for the uniform electron gas is taken 
from Ref. 

resulting E c from Eq. (|21|l is independent of the choice 
of w A e . However, when approximations are made some 
"paths" can give much better results than others [fj. As 
we shall see, this is the case with the present approach. 

We build an Overhauser-like potential for interaction 
w A e (to be added to v^) as 

v^\r 12] r s ) = vi(r 12 )- [ n^ e (|r-r 12 |) dr. (22) 

J\r\<r, 

That is, the average density n of Eq. (fTTjl (and thus 
the average r s ) is kept fixed to mimic the fact that the 
one-electron density does not change along the adiabatic 
connection. The modified interaction u A e is screened by 
a sphere of radius r s and of positive uniform charge of 
density n that attracts the electrons with the same mod- 
ified interaction. This attractive background approxi- 
mates the effect of the external potential V x on /. 

A. Linear adiabatic connection 

If we choose v x e = A/r 12 we simply obtain 
w^' A (ri 2 ;r s ) = \v^(r 12 ;r s ), where (r 12 ;r s ) is 
given by Eq. {2J| . 

The results for (V e x e ) - (V ee ) KS for He and Ne 8 + are 
shown in Fig. 0] and are compared with the "exact" 
ones of Ref. jij. The correlation energy E c can be 
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FIG. 3: The real space analysis of the correlation part of 
the expectation value of V ee : the area under each curve gives 
(V ee ) — (V e e)xs (see also Eq. Q) and Fig. [3Jl. Our results for 
He and Ne 8+ are compared with the "exact" ones and with 
the LDA result (the hole for the uniform electron gas is taken 
from Ref. 0). 



calculated as the area under each curve. We obtain 
E c = -0.052 Hartree for He and E c = -0.053 Hartree 
for Ne 8+ , to be compared with the corresponding 'exact' 
results, -0.042 and -0.045, respectively. 



B. A nonlinear adiabatic connection 

As shown by Figs.nj^and TableU the Overhauser-like 
potential gives accurate results for the short-range part of 
fc(r 12)- We can thus expect to obtain better correlation 
energies from the adiabatic connection formalism if we 
choose a modified interaction i> A e that is able to separate 
long-range and short-range contributions, like the "erf" 
interaction [el \M IH^ 



««(ri 2 ) 



crf(Ari2 
ri2 



With this choice, Eq. I|21f) becomes 

/>oo />oo 

E c [n] = dX dr 12 Anr 2 2 ^{r 12 ) 2 



(23) 



(24) 



For large A, when we are approaching the physical sys- 
tem, the gaussian factor e~ A Vl2 in Eq. I|24(l quenches the 
long-range contribution of / A to the energy integrand. At 
the KS end of the adiabatic connection, when A — > 0, the 
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FIG. 4: Correlation part of {V* e ) along the linear adiabatic 
connection for He and Ne 8+ . Our results are compared with 
the "exact" ones of Ref. |4^|. The area under each curve gives 
the correlation energy E c of standard Kohn-Sham theory. 



interaction, and thus / A , become small, so that the con- 
tribution to E c coming from A- values for which the long- 
range part of / A is not quenched is moderate. Moreover, 
the function / A is correctly normalized to zero so that for 
A — > 0, not only is f£ small, but also the integral itself 
vanishes. In the linear adiabatic connection of Eq. i|12[l . 
instead, the long-range part of / A plays an important 
role in the energy integrand at all A. Indeed, with this 
nonlinear adiabatic connection we obtain E c = —0.0405 
Hartree for He and E c = —0.0413 for Ne 8+ , much closer 
to the "exact" values with respect to the results from the 
linear adiabatic connection. 

The technical details of this calculation are as fol- 
lows. The potential v e ^' X (ri 2 ; r s ) of Eq. j^Sjl can 
be computed analytically, and is reported in Ap- 
pendix[5] We thus obtained, via Eqs. ©-0, dE^/d\ = 
J dr 12 f£(r 12 ) ^7=e~ A Tl2 for 23 values of A between and 

20 for He, and between and 100 for Ne 8+ . We then fit- 
ted our results with the derivative of the following func- 
tional form 



E. 



(1 + 6 2 x 2 Y 



10 



A 

Z° 



(25) 



that has exact asymptotic behaviors for small and large 
A [33. (We have numerical evidence that our results 
fulfill such exact behaviors.) In Fig. we report our 
numerical values for dE^ /d\, together with the derivative 
of the fitting function of Eq. I|25|l . The parameters and 
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FIG. 5: The derivative dE£/d\ along the nonlinear adiabatic 
connection defined by Eqs. I231 - H24|l for He and Ne 8+ . Our 
results are compared with the derivative of the fitting function 
of Eq. I|25jl . The area under each curve from zero to oo gives 
the correlation energy E c of standard Kohn-Sham theory. 



the r.m.s of residuals are reported in Table ITU The KS 
correlation energy E c [n] is then given by 03 /b 10 . 

The accuracy of our results with the "erf" adiabatic 
connection is ofparticular interest for the method of 
Refs. |al2Sl2al23i which combines multideterminantal 
wavefunctions (configuration interaction, CI) with den- 
sity functional theory ("CI+DFT"). In such approach, 
instead of the KS system, one choses a reference system 
of partially interacting particles, usually with the poten- 
tial of Eq. (|23[l . This model system is treated with a 
multideterminantal wavefunction, in a CI fashion, that 
allows to treat near-degeneracy effects. The remaining 
part of the energy is calculated via a density functional, 
that needs to be approximated. The larger A, the larger 
is the energy fraction treated with the CI calculation, 
and thus the larger is the computational cost. The corre- 
lation e nergy fu nctional that needs to be approximated 

is nmnf 



2,y fui 



E c [n] = E c [n]-E x c [7 



(26) 



and can be rewritten as 



— A 



KW= I dX ' I dr 12 4nr 2 12 f* (n 2 ) ^-c 

(27) 

Thus, only the short-range part of f£ contributes to the 

functional £7 c [n], and we expect to get accurate results 
with the present approach. Indeed, this is the CctSG, clS 



He 1.2047 2.3253 2.7788 1.5263 4-10~ J 
Ne 8+ 0.3983 0.4711 0.4026 1.2557 9-10~ 6 

TABLE H: Optimal fit parameters and r.m.s of the residuals 
for the derivative of Eq. J2ol . that parametrizes our results 
for dE* /d\ along the nonlinear adiabatic connection defined 
by Eqs. See also Fig. 03 
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FIG. 6: A nonlinear adiabatic connection that separates long- 
and short-range effects: difference between the correlation en- 
ergy E c of the physical system (with full interaction 1/r) and 
the correlation energy of the system with partial interac- 
tion erf (A r)/r, for He and Ne 8+ . Our results are compared 
with the "exact" ones of Ref. 361 . 



shown in Fig. EI where we compare our results as a func- 
tion of A with the "exact" ones of Ref. [3(j. The error 
is less than 0.5 mHartree for A > l/r s , that is a very 
reasonable choice for the value of A to be used in the 
CI+DFT method of Refs. IE |M IM Ezl . 



VI. CONCLUSIONS AND PERSPECTIVES 

In this work we have started to explore the possibil- 
ity of solving simple radial equations to generate realistic 
spherically- and system-averaged electronic pair densities 
/(^12) for nonuniform systems. With a simple approxi- 
mation for the unknown effective electron-electron inter- 
action that appears in our formalism, we have obtained, 
for two-electron atoms, results that are in fair agreement 
with those coming from accurate variational wavefunc- 
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tions (Figs, dEU and Table P). We have then extended 
our approach along a nonlinear adiabatic connection and 
obtained Kohn-Sham correlation energies whose error is 
less than 4 mHartrees, and short-range-only correlation 
energies whose accuracy is one order of magnitude better 

(Fig.ES). 

In Sec. ^ we have introduced a general formalism for 
many-electron systems that will be further tested in fu- 
ture work. So far we can say that this formalism, com- 
bined with simple physical approximations, works very 
well for two completely different systems: the uniform 
electron gas (23|, |24], [23 and the He series. We think that 
this fact makes the method promising. 

To fully develop the approach described in this paper, 
many steps have to be performed. First of all, the KS 
part of the effective e-e potential, (r^) of Eq. (jTHjl . 
should be also approximated, to make the extension to 
many-electron systems practical. The correlation part of 
the effective e-e potential can be improved, in analogy 
with the recent developments for the uniform electron- 
gas case [3 H3 • It should then be possible to construct 
a self-consistent scheme (OEP-like) that combines the 
Kohn-Sham equations with the correlation energy func- 
tional arising from our approach [Eqs. at different 
coupling strengths A, plus Eq. <|T2*)l or Eq. JUJl]. With 
respect to traditional DFT calculations, this combined 
scheme would have the advantage of yielding not only 
the ground-state one-electron density n(r) and energy E, 
but also the spherically- and system-averaged pair den- 
sity f(r '12), thus allowing to calculate expectation values 



of two-body operators that only depend on the electron- 
electron distance. The combination of our approach with 
the CI+DFT method of Refs. HEHUll could be also 
implemented and, in view of the results of Fig. it is 
even more promising. We are presently working in all of 
these main directions. 
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APPENDIX A: OVERHAUSER-LIKE 
POTENTIAL FOR THE ERF INTERACTION 

The evaluation of Eq. (122[l with the interaction v£ e — 
erf(Ari 2 )/ri2 gives 



Ov,\, _ \ 



u(s,fj,) 



(Al) 



where s = ryij^s, M = ^ r s, and 



J 



u(s,fj,) = 



erf (/x s) 



1 



8 \pii s /j, 3 



l+{-2 + S + s 2 ) ^ + e- i ^ 2 (-l + (2 + s-s 2 ) /i 2 ) 



V^fj, 3s + 2{l- S y (2 + s) n 2 erf \p (l-s)]+^/j, -3 s + 2 (2 - s) (1 + s) V erf [/i (1 + s)} }. (A2) 



[1] W. Kohn, Rev. Mod. Phys. 71, 1253 (1999). 

[2] A.E. Mattsson, Science 298, 759 (2002). 

[3] C. Fiolhais, F. Nogueira, and M. Marques (eds.), A 

Primer in Density Functional Theory (Springer- Verlag, 

Berlin, 2003). 

[4] W. Kohn and L.J. Sham, Phys. Rev. 140, A1133 (1965). 

[5] J. P. Perdew and K. Schmidt, in Density Functional The- 
ory and Its Applications to Materials, edited by V. Van- 
Doren et al. (AIP, NY, 2001), and references therein. 

[6] A. Savin, F. Colonna, and R. Pollet, Int. J. Quantum 
Chem. 93, 166 (2003), and references therein. 

[7] V. Sahni, Quanta! Density Functional Theory (Springer- 
Verlag, Berlin, 2004). 

[8] J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 
Lett. 77, 3865 (1996); ibid. 78, 1396 (1997); A.D. Becke, 



Phys. Rev. A 38, 3098 (1988); J. Chem. Phys. 84, 4524 
(1986); C. Lee, W. Yang, and R.G. Parr, Phys. Rev. B 
37, 785 (1988); J.P. Perdew, J.A. Chevary, S.H. Vosko, 
K.A. Jackson, M.R. Pederson, D.J. Singh, and C. Fiol- 
hais, Phys. Rev. B. 46, 6671 (1992); ibid. 48, 4978 (1993). 
[9] J. Tao, J.P. Perdew, V.N. Staroverov, and G.E. Scuseria, 
Phys. Rev. Lett. 91, 146401 (2003). 
[10] J.P. Perdew, A. Ruzsinszky, J. Tao, V.N. Staroverov, 
G.E. Scuseria, and G.I. Csonka, J. Chem. Phys., to ap- 
pear. 

[11] M. Seidl, J.P. Perdew, and S. Kurth, Phys. Rev. Lett. 
84, 5070 (2000). 

[12] see, e.g., S. Ktimmel and J.P. Perdew, Phys. Rev. B 68, 
035103 (2003); W. Yang and Q. Wu, Phys. Rev. Lett. 89, 
143002 (2002); R. J. Magyar, A. Fleszar, and E. K. U. 



9 



Gross, Phys. Rev. B 69, 045111 (2004); M. Griming, O. 
V. Gritsenko, and E. J. Baerends, J. Chem. Phys. 118, 
7183 (2003). 

[13] C.A. Coulson and A.H. Neilson, Proc. Phys. Soc. London 
78, 831 (1961). 

[14] J. Cioslowski, B.B. Stefanov, A. Tan, and C.J. Umrigar, 

J. Chem. Phys. 103, 6093 (1995). 
[15] J. Cioslowski and G. Liu, J. Chem. Phys. 109, 8225 

(1998). 

[16] E. Valderrama, J.M. Ugalde, and R.J. Boyd, in Many- 
electron densities and reduced density matrices, edited 
by J. Cioslowski (Kluwer Academic/Plenum Publishers, 
New York, 2000). 

[17] E.R. Davidson, Reduced Density Matrices in Quantum 
Chemistry (Academic Press, New York, 1976). 

[18] A.J. Coleman and V.I. Yukalov, Reduced Density Ma- 
trices: Coulson's Challenge (Springer- Verlag, New York, 
2000). 

[19] see, e.g., G. Mazzone, F. Sacchetti, and V. Contini, Phys. 

Rev. B 28, 1772 (1983); C. Petrillo and F. Sacchetti, ibid. 

51, 4755 (1995). 
[20] A. K. Rajagopal, J. C. Kimball, and M. Banerjee, Phys. 

Rev. B 18, 2339 (1978); X.-Y. Pan and V. Sahni, J. 

Chem. Phys. 119, 7083 (2003). 
[21] The difference in kinetic energy can be obtained from 

the adiabatic connection formula that can be viewed as 

a generalization of the above argument, and is treated in 

Sec. |V| 

[22] A.W. Overhauser, Can. J. Phys. 73, 683 (1995). 
[23] P. Gori-Giorgi and J.P Perdew, Phys. Rev. B 64, 155102 
(2001). 

[24] B. Davoudi, M. Polini, R. Asgari, and M.P. Tosi, Phys. 

Rev. B 66, 075110 (2002). 
[25] M. Corona, P. Gori-Giorgi, and J.P. Perdew, Phys. Rev. 

B 69, 045108 (2004); I. Nagy, R. Diez Muino, J.I. Juar- 

isti, and P.M. Echenique, Phys. Rev. B 69, 233105 

(2004). 

[26] P. Ziesche, Phys. Lett. A 195, 213 (1994); M. Levy and 

P. Ziesche, J. Chem. Phys. 115, 9110 (2001). 
[27] A. Gonis, T.G. Schulthess, J. van Ek, and P.E.A. Turchi, 

Phys. Rev. Lett. 77, 2981 (1996); A. Gonis, T.G. 

Schulthess, P.E.A. Turchi, and J. van Ek, Phys. Rev. 

B 56, 9335 (1997). 
[28] A. Nagy, Phys. Rev. A 66, 022505 (2002). 
[29] F. Furche, Phys. Rev. A 70, 022514 (2004). 
[30] J.M. Soler, Phys. Rev. B 69, 195101 (2004). 
[31] M. Levy, Proc. Natl. Acad. Sci. U.S.A. 76, 6062 (1979). 
[32] E. Lieb, Int. J. Quantum Chem. 24, 243 (1983). 
[33] J. Harris and R. Jones, J. Phys. F 4, 1170 (1974); D.C. 



Langreth and J.P. Perdew, Solid State Commun. 17, 
1425 (1975); O. Gunnarsson and B.I. Lundqvist, Phys. 
Rev. B 13, 4274 (1976). 

[34] W. Yang, J. Chem. Phys. 109, 10107 (1998). 

[35] A. Savin, in Recent Developments and Applications of 
Modern Density Functional Theory, edited by J.M. Sem- 
inario (Elsevier, Amsterdam, 1996); T. Leininger, H. 
Stoll, H.-J. Werner, and A. Savin, Chem. Phys. Lett. 
275, 151 (1997); R. Pollet, A. Savin, T. Leininger, and 
H. Stoll, J. Chem. Phys. 116, 1250 (2002). 

[36] R. Pollet, F. Colonna, T. Leininger, H. Stoll, H.-J. 
Werner, and A. Savin, Int. J. Quantum Chem. 91, 84 
(2003); J. Toulouse, private communication. 

[37] J. Toulouse, F. Colonna, and A. Savin, Phys. Rev. A 70, 
062505 (2004). 

[38] B. Davoudi, R. Asgari, M. Polini, and M.P. Tosi, Phys. 
Rev. B 68, 155112 (2003); R. Asgari, B. Davoudi, and 
M.P. Tosi, Solid State Comm. 131, 301 (2004). 

[39] R. Asgari, B. Davoudi, M. Polini, and M.P. Tosi, in 
preparation. 

[40] D.E. Freund, B.D. Huxtable, and J.D. Morgan III, Phys. 
Rev. A 29, 980 (1984). We used an improved version 
(provided to us by C. Umrigar) of the accurate varia- 
tional wavefunctions described in this work to obtain one- 
electron densities n(r) and functions /(ri2). See also C.J. 
Umrigar and X. Gonze, Phys. Rev. A 50, 3827 (1994), 
and Ref. [3. 

[41] P. Ziesche, Phys. Rev. B 67, 233102 (2003); P. Ziesche, 
K. Pernal, and F. Tasnadi, Phys. Status Solidi B 239, 185 
(2003). In these papers interesting exact properties of the 
"Overhauser geminals" are derived, using the interacting 
momentum distribution to define the geminal occupancy. 
Notice however that the uniform-electron gas equivalent 
of our Eqs. JSJ-Q employ the non-interacting momen- 
tum distribution, as it has been done in Refs. |23l l24ll25| . 

[42] F. Colonna and A. Savin, J. Chem. Phys. 110, 2828 
(1999). 

[43] Q. Zhao, R. C. Morrison, and R. G. Parr, Phys. Rev. 

A 50, 2138 (1994); R. van Leeuwen and E. J. Baerends 

Phys. Rev. A 49, 2421 (1994). 
[44] J.P. Perdew, A. Savin, and K. Burke, Phys. Rev. A, 51, 

4531 (1995); K. Burke, J.P. Perdew, and M. Ernzerhof, J. 

Chem. Phys. 109, 3760 (1998); E. Valderrama and J.M. 

Ugalde, Int. J. Quantum Chem. 86, 40 (2002). 
[45] X. Fradera, M. Duran, E. Valderrama, and J.M. Ugalde, 

Phys. Rev. A 62, 034502 (2000). 
[46] P. Gori-Giorgi and J.P. Perdew, Phys. Rev. B 66, 165118 

(2002). 



